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As waves propagate into shallower waters, their profile changes fi'om a Gaussian 
distribution to an increasingly non-Gaussian profile. The non-Gaussian properties of 
waves in finite water depths are analyzed through the use of a probability density fimction 
defined in closed form and that can be solved by the application of wave displacement 
data. Initially, the applicability of the probability density function is verified using wave 
data records obtained fi'om the Coastal Engineering Research Center (CERC) at Duck, 
North Carolina. A broad range of sea conditions and water depths are represented to 
ensure the probability density function correlates well with actual histogram data over 
varied conditions. With the applicabihty verified, the probability density function is used 
to define criteria for the boimdary where the wave field can no longer be considered 
Gaussian. The ability to determine when a shoaling wave profile can no longer be 
considered Gaussian is of paramount importance in order to ensure proper wave theories 
are being applied. This boimdary between Gaussian and non-Gaussian wave profiles will 
be determined in terms of significant wave height and water depth. 
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Presently, wave data records are required to define the probability density function 
used in this study. Trend analyses were conducted on the parameters of the probability 
density function to determine if sinq)lifications can be made to its form and to see if it is 
possible to define the probabUity density function in terms of the local conditions, sea 
severity and water depth, without the need for wave data records. 
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CHAPTER 1 
INTRODUCTION 



As wind generated waves transition from deep water to shallow water, a significant 
change is evident in the wave profile. As an example. Figure 1-1 shows portions of wave 
records obtained by the Coastal Engineering Research Center, (CERC), during the 
ARSLOE Project. Figure l-l(a) is the wave profile recorded at a distance of 60 meters 
ofl&hore in a shallow intermediate water depth condition. Figure 1-1 (b) is the wave 
profile taken at the same time but at a distance of 12 km offshore in a deep water 
condition. Deep water wave profiles in the sea characteristically reflect a Gaussian (linear) 
random process. The height of the wave crests are probabUistically equal to the height of 
the wave troughs. As waves move from deep to shallow water, the wave profile 
transforms from a Gaussian random process to a non-Gaussian (nonlinear) random 
process in which the wave profile shows a definite excess of high crests and shallow 
troughs. 









a) Gage 615, Severe Sea Condition b) Gage 625, Mild Sea Condition 



Figure 1.1 Wave Profile Examples 
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This remarkable transformation in the wave profile, as an irregular sea transitions 
from deep to shallow water, has also been evidenced in laboratory experiments (Doering 
and Donelan 1993). The transformation of wave profiles can be attributed to the non- 
linear wave-wave interactions augmented primarily by decreasing water depth. The 
“degree” of non-linearity a wave field has at a particular site of interest in the coastal 
region appears to depend largely on sea severity and water depth. 

Although there has been awareness that CJaussian-based solutions do not always 
represent accurate results in the coastal region, clarification of wave characteristics from 
deep to shallow water has been a subject of considerable interest in ocean and coastal 
engineering. There have been few studies which have offered non-Gaussian solutions 
describing the probability density fimction of waves in coastal regions. A commonly 
known study was by Longuet-Higgins (1963), which statistically predicts coastal wave 
properties using a power series ejqjansion. This probability distribution is based essentially 
on the Gram-Charlier series distribution. A probability density fimction developed by Ochi 
and Ahn (1994a) appears to be applicable throughout the coastal region from deep to 



shallow water and is defined in closed form Equation (1). 



f(x) = 



1 






exp 



1 



2(yaa,y 



-(] - ya//. - exp(-yax))^ - yax 



(I) 



where, y = { 



L28 for.x > 0 



3.00 for: x < 0 

This is the probability density fimction that wiU be used as the basis to analyze 
coastal wave profiles for this study. 



1.1 Purpose of the Study 



The overall objective of this study involves two components. The first component 
will be to vabdate the applicability of the probability density fimction developed for waves 
in finite water depth (Ochi and Ahn 1994a) and to analyze the transition point where wave 
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profiles can no longer be considered Gaussian. In order to accomplish the first objective of 
this study, a statistical analysis was carried out on wave time histories obtained dming a 
storm by the Coastal Engineering Research Center (CERC) during the Atlantic Ocean 
Remote Sensing Land-Ocean Experiment (ARSLOE) Project in 1980. To interpret the 
transformation of waves fi^om a Gaussian to a non-Gaussian random process, wave 
records at seven locations were taken where the water depths ranged fi'om two meters to 
twenty-one meters. In total, approximately 1000 wave records were analyzed with each 
wave record consisting of 4800 water displacement measurements. The individual wave 

I 

displacement measurements were recorded at a fi'equency of 4 Hertz in twenty minute 
intervals. 

Upon confirmation of the probabihty density fimctions’ apphcabihty, the histogram 
data will be compared closely with the corresponding normal distribution to determine the 
conditions under which the wave profile can no longer be considered Gaussian. In 
particular, general relationships to sea severity and water depth will be determined. A 
criteria wih be estabhshed which will define the limit to when the Gaussian wave profile 
becomes non-Gaussian. This limit will be based on the nonlinearity parameters analyzed 
fi'om the wave profile data. The second objective of the study will then be to analyze the 
parameters of the probability density function. The three parameters, a, p* and a* wiU be 
statistically analyzed in both dimensional and non-dimensional form Trends based on the 
analysis wiU be described and relationships linking these parameters of the probability 
density function to site conditions wiU be developed. 

This study consists of six chapters. Chapter 2 provides a general re’view on the 
development of probability density functions apphcable to coastal regions. Chapter 3 
describes the ARSLOE Project of 1980 and addresses the statistical approach, taken in this 
study. The data used as the basis for this study, were obtained fiom the ARSLOE Project. 
This chapter will also describe the local field conditions at CERC during the time the data 
were collected. Chapter 4 addresses the probability density function applicability and 
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analysis. This chapter has two distinct sections. The first section reviews the finite water 
depth probabihty density fiinction (Ochi and Ahn 1994b) and verifies the applicabihty of 
the probabihty density fiinction. The second section discusses sea severity and water 
depth and the transition between Gaussian and non-Gaussian wave profiles. Chapter 5 
relates the non-dimensionalized parameters of the probabihty density fiinction, ac, p./a 
and a./a, to sea severity, (a), and water depth in both dimensional and non-dimensional 
form. The parameters will also be compared to each other and relationships determined. 
This chapter describes trends evident based on statistical analysis. These trends are 
quahtative and are intended to give a better xmderstanding of the interaction the wave 
profile has with its environment as it propagates in shaUow waters. FinaUy, Chapter 6 
provides the summary and conclusion for this study. 



CHAPTER! 
LITERATURE REVIEW 



As stated in the Introduction, the profile of waves in finite water depths is quite 
different fi'om that observed in deep water. From the stochastic view point, waves in finite 
water depth are categorized as a non-Gaussian random process for which a precise 
description of the process does not exist. This is in contrast to the Gaussian random 
process which can be defined by a rigorous mathematical formulation. Although 
considerable attention has been given to Gaussian waves in deep water, little information 
is available on the non-Gaussian characteristics of shallow water waves. 

One of the most commonly known presentations describing non-Gaussian waves is 
the Gram-CharUer series probability density fimction. This density function is given in a 
series consisting of an infinite number of terms. Longuet-Higgins (1963) theoreticaUy 
showed that the statistical properties of nonhnear waves can be presented by the Gram- 
Charher series distribution. Bitner (1980) found reasonable agreement in a comparison 
between the Gram-CharUer series distribution and the histogram constructed fi'om wave 
data obtained in a relatively shallow water area. 

On the other hand, Ochi and Wang (1984) analyzed more than 500 wave records 
obtained by the Coastal Engineering Research Center (CERC) at Duck, North Carolina 
during the ARSLOE Project (1980). They found that while the statistical properties of 
waves at any water depth can be well represented by the Gram-Charher series distribution, 
there is a drawback to applying this density function in practice. The probability density 
function becomes negative for large negative wave displacements and this negative portion 
does not vanish with increasing number of terms in the series distribution. Because of this 
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difficulty, the transition property of wave characteristics from gaussian to non-gaussian 
distribution can not be well clarified. They derived, however, a limiting sea severity for a 
specified water depth below which wind-generated coastal waves may be considered to be 
a Gaussian random process. 

Another approach for describing the surface elevation of non-Gaussian random 
waves is to assume that individual waves can be expressed as a Stokes wave. There 
appears to be some reservation to applying this approach for representing the statistical 
distribution of the wave profile. This method presupposes a preliminary form of the wave 

t 

profile before estabUshing its probabihty density. This supposition may impose 
unfavorable bias on the outcome. 

Tayfim (1980) derived the probabihty density fimction based on the Stokes profile 
approximation up to the second order. Huang et al. (1983) developed the probabihty 
density fimction considering ah terms up to the third order. These probabihty density 
fimctions are superior to the Gram-Charher series distribution in that they do not have any 
negative portion over the entire range of wave displacement. Comparison of these 
probabihty density functions and wave data show good agreement. 

Langley ( 1987) apphed a probabihty density function representing the response of 
a non-linear system developed by Kac and Siegert (1947), to nonlinear random waves. 

This probabihty density function is given in the form of a series and can only be evaluated 
numericaUy. Ochi and Ahn (1994a) modified the Kac Siegert density function so that it is 
presented in closed form with parameters determined through spectral analysis. They 
apphed the probabihty density function to non-Gaussian coastal waves. In the apphcation 
process, they decomposed the spectra into linear and nonlinear components to clarify the 
degree of non-linearity involved in shaUow water waves (Ochi and Ahn 1994b). The 
probabihty density function does not have any negative portion in its entire domain and 
appears to represent the statistical distribution of the wave profile in any sea state during a 
storm in deep, intermediate and shaUow water. 
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Studies concerned with the non-Gaussian properties of waves in finite water depth 
include Thornton and Guza ( 1983) who claimed that the results of their analysis of field 
data obtained at Torrey Pines Beach, California, show that the Rayleigh probability 
distribution can be apphed to represent wave heights observed fi'om shallow to deep (10.8 
meter) water. This inches that all observed waves can be considered a Gaussian random 
process. However, the sea severity in relation to water depth at all locations in the 
nearshore zone \\iiere their measurements were taken was too mild and below the limit for 
waves to be a non-Gaussian random process. This subject will be discussed in detail in 
this present study. 

In closing this chapter, it is important to note that there still remains a significant 
amount of work to do in refining our abihty to accurately describe irregular seas in finite 
water depths and over a wide range of sea severity. 



CHAPTERS 

BACKGROUND / STUDY APPROACH 
3. 1 Wave Record Data 

An experiment jointly hosted by the U.S. Army Coastal Engineering Research 
Center (CERC) and the National Ocean Survey (NOS) was conducted from 06 October 
1980 to 30 November 1980. An aspect of this Atlantic Ocean Remote Sensing Land- 
Ocean Experiment (ARSLOE) involved ocean wave studies implemented near the 
CERC’s Field Research Facility (FRF) at Duck, North Carolina. In situ wave measuring 
devices were deployed in mean water depths ranging from 1.35 meters to 24.4 meters and 
were placed in a linear array extending from 60 meters offshore to 12,000 meters offshore. 
Figure 3.1. 

3.1.1 Data Acquisition 

The wave displacement data analyzed in this study came from a total of eight wave 
measuring devices. The FRF has a 550 meter research pier oriented perpendicular to the 
local shore which was equipped with seven Baylor resistance-type wave staffs. These 
Baylor wave staffs function by determining the change in electrical resistance experienced 
by changes in water elevations. Round-trip travel time of an electrical current from an 
above-water emitter to the water surface is measured. Following along the line of the 
pier off-shore was a Waverider Buoy Gauge. This is a surface-following class buoy. It 
functions by double-integrating the output of an accelerometer held vertical by a pendulum 
to produce the time history of buoy heave motion. The statistics of the buoy heave then 
represent the wave displacements. In the Waverider buoy, the pendulum is contained 
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entirely wi thin the spherical float so the buoy does not have to be stabilized. These 
devices recorded both wave displacements and wave periods but only the displacement 
data were analyzed for this study. The wave gauges recorded water surface fluctuations 
simultaneously at all eight locations at a fi'equency of 4 Hertz over twenty minute 
intervals. This resulted in 4800 individual wave displacements recorded for each twenty 
minute segments continuously. Due to the signiScant amount of data obtained during this 
unique experiment, only a sub-set of the data was analyzed. This chosen data sub-set 
began at 1200 on 23 October 1980 and ended at 1315 on 25 October 1980. This two day 
segment was unique in that a severe storm occurred fi^om 23-25 October which was the 
most significant meteorological event recorded diuing the experiment. Wind speeds on 
the order of fifteen meters per second were recorded during the severe storm. This two 
day sub-set of data was used as the basis for this study since it included the entire range of 
wave profiles experienced at each of the eight wave gauge locations. From this two day 
sub-set, a total of nearly 5 million individual wave displacements were analyzed and the 
wave profile data covered the entire range of sea conditions present during the complete 
56 day ARSLOE Project. 

3.1.2 Wave Profile Transformations 



The data collected at the FRF during ARSLOE is unique and impressive. It is very 
well suited for analyzing changes in the wave field as the waves propagate through 
shoaling depths. The bathymetry at the experiment location is relatively simple showing 
straight and parallel contours. This provided favorable conditions which maximized 
control over the environmental impact on the quaUty of the data received. The wave 
profile data provide a matrix of wave profiles allowing both spatial and temporal analyses 
to be performed. 
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Figure 3.1 Coastal Engineering Research Center Field Research Facility 
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A wide range of wave profiles was observed during the e?q)eriment. Figures 3.2- 
3.8 represent a nuld sea condition recorded on 10/23 @2155 for all finite water depth 
gage locations shown in Figure 3.1. Figures 3.9-3.15 represent a severe sea condition 
recorded at the height of a storm on 10/25 @ 1315 for the same gage locations. It is 
evident from the sample wave data that the wave field transitions as h shoals from deep to 
shallow water. The wave data representing the mild sea condition at Gage 625, Figure 
3.8, are characterized by an almost symmetrical distribution of positive and negative wave 
profile displacements. This is characteristic of Gaussian random processes. As the wave 
shoals into shallower water, the wave profile transforms showing progressively higher 
peaks and shallower troughs. This process is clearly non-Gaussian and must be analyzed 
differently. It is also clear from observing these wave records that the degree of non- 
linearity increases as the relative water depth decreases and as sea severity increases. This 
Lmphes that analysis of this nonlinearity must include a dependence on water depth as it 
relates to the magnitude of the waves’ displacement (sea severity). 



12 




Figure 3-2 Gage 615 - 10/23 @ 2155 (d=1.74m) 






Figure 3-3 Gage 635 - 10/23 @ 2155 (d=2.83m) 






Figure 3-4 Gage 645 - 10/23 @ 2155 (d=4.31m) 






Figure 3-5 Gage 655 - 10/23 @ 2155 (d=6.0m) 
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Figure 3-6 Gage 665 - 10/23 @ 2155 (d=6.08m) 



Figure 3-7 Gage 675 - 10/23 @2155 (d=6.82ni) 







Figure 3-8 Gage 625 - 10/23 @ 2155 (d=9.66m) 
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Figure 3-9 Gage 615 - 10/25 @ 1315 (d=2.23m) 



Figure 3-10 Gage 635 - 10/25 @ 1315 (d=3.32m) 



Figure 3-11 Gage 645 - 10/25 @ 1315 (d=4.66m) 







Figure 3-12 Gage 655 - 10/25 @ 1315 (d=6.44m) 



Figure 3-13 Gage 665 - 10/25 @ 1315 (d=6.52m) 
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Figure 3-15 Gage 625 - 10/25 @ 1315 (d=10.04m) 
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3.2 Nearshore Zone Characteristics 

All wave fields propagating in finite water depths are influenced by the attributes 
of the local sea environment. These attributes include water depth and beach slope. An 
awareness of the characteristics of these attributes is inq)ortant to understand the type of 
influence they have on the wave field. 

3.2.1 Discussion of Water Depth Redons During ARSLOE 

Based on linear theory, there are three water depth regions which are defined using 

the velocity potential for a progressive wave. The velocity potential is defined as 

HgT co^{k{d z)) 2n*t 

The distinct water depth regions are determined by the asymptotes of the hyperbolic 
function. Deep water wave data obtained at gage location 710 and shown in Table 3. 1 
allow these regions to be defined for the mild sea condition and the severe sea condition. 

In general, the deep water wave condition is defined as kd > n, where k is the wave 
number and h is the water depth. This condition occurs at d = 2 Im for the mild sea 
condition on 10/23 @2155 and d = 43m for the severe sea condition on 10/25 @ 1315. 
The shallow water condition is defined as kd < tc/ 10. This condition occurs at d = 0.6m for 
the mild sea condition and d = 1.3m for the severe sea condition on 25 October. The 
intermediate water condition occurs between these boimds. All wave profile data 
obtained from ARSLOE are located in the intermediate water depth region. 



Table 3.1 Deep Water Wave Statistics for Mild and Severe Sea Conditions 



Date/Time 


Zero Crossing Period, T 


Deep Water Wavelength, X.n 


10/23 @ 2155 


5.17 sec. 


41.715 m 


10/25 @ 1315 


7.43 sec. 


86.157 m 
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3.2.2 Discussion of the Surf Zone During ARSLOE 

An integral aspect of any study ^\1lich covers the nearshore region must analyze 
and consider the surf zone. The surf zone is a transient entity which varies in width off- 
shore depending on the local sea conditions. These local conditions include the beach 
profile, sea severity and direction of wave field. In order to define the surf zone in this 
study, wave field data were analyzed which represented the range of sea conditions 
experienced over the two day period. The first condition, representing the mildest sea 
condition, was taken as 2155, 23 October 1980. The second condition, representing the 
most severe sea condition was taken as 1315, 25 October 1980. 

Linear wave theory will be used to determine the effects of shoaling as the wave field 
enters shallow water. Refraction will not be considered since this process tends to reduce 
the shoaling wave heights and thus move the surf zone shoreward. This is shown by the 
definition of the refraction coeflBcient, where Go and 0 are the deep water 

V cos^ 

wave angle and the finite water depth wave angle respectively. This ratio will always be 
less than one causing a reduction in the resulting wave height. In order to ensure that 
conservative results are obtained, refraction will be neglected in this study. The beach 
profile at the ARSLOE project site shows shghtly varying beach slopes so a wave 
breaking criteria which considers the beach slope will be used. Weggel (1972) 
reinterpreted laboratory results and showed a dependency of breaker height on beach 
slope. The following equations were used to define the surf zone locations: 



Conservation of Energy: H ^ = H 




HI 
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Breaker Criterion; , where k is defined by Weggel as follows: 

K = b(m) - with a(m) = 43.8(L0 - 



gT 



6(m)=L56(L0 + g 



As the beach slope, m approaches zero, this equation reduces to k = 0.78. 



The Group Velocity, Cg = 



2kd 



sinh(2^ri) 



will be determined based on the 



dispersion relation. 



Itt 

~f 



= gk tanh(AzZ) , where L = wave length, T = wave period 



(assumed constant), g = gravitational acceleration, k = the wave number and d = local 
water depth. 

In order to standardize the location of the surf zone, the root-mean-square wave 

n ^ r 

height, Hms , will be used to represent the irregular wave field. H = — V ^ 

V ^ .=1 



where H; is defined as the wave height of the i* wave of the irregular wave field. The surf 
zone which results firom the use of H„ns will represent the approximate “center of mass” 
location where the most significant set of individual waves in the wave field begin to 
break. The surf zone calculations were based on a root-mean-square wave height that was 
determined fi-om the deep water wave displacement data. This single wave height was 
used to represent the irregular sea and shoaled landward. Since the real sea will have a 
significant number of waves with heights larger and smaller than the root-mean-square, the 
actual breaking zone is an obscure entity. This technique for determining the surf zone 
locations is adequate to make general comparisons when analyzing the spatial changes in 
the nonlinearity of the irregular sea in finite water depths. The surf zone locations during 
the mildest and most severe sea conditions over the two day period are shown in Table 
3.2. 



In Table 3.2 it is shown that during the mild sea condition, the breakpoint, indicating 
the beginning of the surf zone, is located 70 meters from the shoreline in 1.5 1 meters of 
water. This breakpoint is the location where a wave having a height equivalent to the 
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root-mean-square wave height of the wave field will break. During the severe sea 
condition, this breakpoint occurs 1 10 meters fi'om the shoreline in 3.32 meters of water. 
Realistically, waves would be observed breaking considerably seaward of the breakpoints 
determined in the respective sea conditions due to the spectrum of wave heights present in 
an irregular sea. The root-mean- square wave height breakpoint is used in order to 
standardize the position of the surf zone for camparison purposes. 



Table 3.2 Surf Zone Locations for Mild and Severe Sea Conditions 



Sea Condition 


Depth where a wave v»4th 


Start of Surf Zone 




height Hnns 7 Breaks 


-Dist. from Shore- 


Mild - 10/23 @ 2155 


1.51 m 


70 m 


Severe - 10/25 1315 


3.32 m 


110m 



3.3 Non-Gaussian Probabihtv Density Function 

The probabihstic approach used in this study wUl make use of the probability 
density fimction appHcable to non-gaussian random processes developed by Ochi and Ahn 
(1994a). This probabihty density fimction appears to represent non-Gaussian waves in 
shallow and intermediate water depths as well as Gaussian waves in deep water. 

3.3.1 Non-Gaussian Probabihtv Density Fimction Backgroimd 

The probabihty density fimction has a backgroimd of the Kac-Siegert solution 
developed for the response of non-hnear systems. The density fimction, y(t), is given as 
the sum of the standardized normal variates and its squares as follows 
y(t)= (p ,2^ + ^jZ/) 

J=1 
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The parameters Pj and X,j are determined by finding the eigenfimction and eigenvalues 
of the integral equation given by 



otherwise, v|/j(co) = 0 

H(coi , ©2 ) = second order frequency response fimction 
S(©) = output spectral density fimction 

As can be seen from the form of this solution, the probability density fimction of y(t), 
can not be foimd without knowledge of the second order frequency response fimction, 
H(©i , ©2 ) and there is no way to evaluate this response fimction for irregular seas. Ochi 
and Ahn showed however, that Kac-Siegerts’ solution can be obtained from the time 
history of the wave record without any knowledge of the second order frequency response 
fimction, H(©i , ©2 ). The probability density fimction thus obtained for describing non- 
gaussian random processes is 




where K(©i , ©2 ) = H(©i , ©2 ) y]S{cOi)S{o}2 ) 

v|/j(©) = orthogonal eigenfimction = { J ¥ = 1 for j = k 





, L28/orx>0 




The parameters, a, and are determined from the cumulants evaluated from the wave 



data. These cumulants are related to the Pj and A,j values of the Kac-Siegert solution by 



the following relationships 



k,= 2: A , 



2 N 



2 N 



^2- 'Yj J ^ J 




2 ^ 



k3 = 6X/0^-A,+8X^^-,etc. 
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In the present study, kj represents the mean water level and is therefore defined by, 
ki = 0. k 2 and ka are equal to the data sample moments, E[x^] and E[x^] respectively. 
Based on this relationship, the parameters, a, p* and can be determined through the 
iterative solution of three equations along with the analysis of a data record for the non- 
Gaussian random process being analyzed. The three equations are 
a + a + p» = E[x] = 0 
CT^. - 2a^a‘*. = E[x^] 

2a a\(3 - = E[x^] 



3.3.2 Probability Density Function Application to Wave Data 

In order to determine values for the three parameters, a, p* and a'* , which are 
used to define the degree of non-linearity a shoaling wave field has at a given location and 
sea state, data records must be analyzed. With the applicability of the probability density 
fimction verified based on the wave data collected dming ARSLOE, the approximately 
1000 data records fi'om 1200 on 23 October to 1315 on 25 October will be analyzed to 
determine values for the three parameters defining the probability density fimction. The 
procedure used during this study to perform this task is as follows: 

- Compute the standard deviation of each twenty minute wave record consisting of 
4800 discrete wave displacement data points E(x ^ ) = ~ • 

Hi 

(E[x] does not need to be computed since it is defined as zero and the data are 
already standardized). 



- Compute the third moment of the same twenty minute wave record 
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- Define the standard deviation of the data record as cr = ■\jE(x ^ ) . 

- Evaluate the three parameters, a, p* and o^. , through iterative solution of the 

foUowing three equations; 

a + a + p. = E[x] = 0 
a^-2aV. = E[x^] 

2a a\(3 - 8a V.) = E[x^]. 

It is important to note that the standard deviation of the data record is used in the 
analysis as an indicator of the sea severity and not the significant wave height. This is 
because the definition of significant wave height is not apphcable to non-Gaussian waves. 

All available data sets covering the seven finite water depth gage locations 
between 1200, 23 October 1980 and 1400, 25 October 1980 were analyzed and a, p» and 
were determined. 

Based on a plot of the individual parameters in the space domain, it was 
determined that a smoothing technique should be used to refine the “a” parameter values. 
This smoothing was done in order to remove the sporadic tertiary effects while ensuring 
the smoothed “a” value does not lose any authenticity. An exponential smoothing 
technique was used which resulted in the lowest mean absolute deviation (MAD) between 
the raw “a” parameter value and the smoothed “a” parameter. The remaining two 
parameters, p. and did not show the same susceptibihty to tertiary effects so no 
smoothing was done. Comparisons between the actual and smoothed “a” parameters are 
shown in Appendix A for all Gage positions. 

3.3.3 Dimensional Analysis of PDF Parameters 

For each of the approximately one-hundred-fifty data sets per gage position, the 
probabUity density fimction parameters, a, p, , a» were determined iteratively. Each of 
the parameters contain discrete dimensions. Based on dimensional analysis, these 
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dimensions are as follows: It is known that k 2 , representing the data sets’ variance has 
dimensions of [L^]. It is also known that ks, representing the data sets’ third moment has 
dimensions of [L^]. It follows from the knowledge of these dimensions that the 
parameters, a, p, , a* must have dimensions [1/L], [L] and [L] respectively. 

In order to minimiz e any dependence on external conditions and to ensure the 
solutions are imiversally applicable, these parameters, a, p. , a* must be made non- 
dimensional. The non-dimensional quantities that will be used in this study are, aa, p./o, 
(5 Jo and a/d. c is defined as the standard deviation of that wave displacement data and d 
is defined as the water depth. 



CHAPTER 4 

PROBABILITY DENSITY FUNCTION ANALYSIS 
4.1 Probability Density Function Applicabilitv Verification 

As stated in the introduction, one of the objectives of this study is to determine 
whether or not the probabihty density function developed by Ochi and Ahn is vahd for 
representing waves in finite water depths. To verify the applicability of this probability 
density function, wave displacement data sets obtained from the ARSLOE experiment 
were used in comparing the histogram of the data record to the probability density 
function. Each data set contained 4800 wave displacement readings taken over a twenty 
minute period. To ensure a thorough verification of the probability density fimction, the 
data records used for the analysis were from all seven “finite” water depth gage locations 
and analysis was done over a wide range of sea state conditions. One analysis was done 
on data from 23 October 1980 at 2155. This represented a mild sea state. The other 
analysis was done on data from 25 October 1980 at 1315. This data set was at the peak 
of the most severe storm experienced during the entire ARSLOE study. 

Figures 4. 1-4. 14 show comparisons between the histogram data and the non- 
Gaussian probability density function. Included also in the figures is the Gaussian 
(normal) probability density function for comparison. Figures 4. 1-4.7 pertain to the mild 
sea condition while Figures 4.8-4. 14 represent the severe sea condition. As can be seen 
from these figures, the probability density function agrees very well with the histogram 
generated from the data record in all cases. This imphes that the probability density 
function derived by Ochi and Ahn can be used with confidence when analyzing non- 
Gaussian random processes such as wave displacements in shallow and intermediate 
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waters. It can also be seen that the probability density iunction approximates the normal 
distribution when the wave field is in the deeper intermediate water condition under mild 
sea conditions. 




Table 4.1 Gage 615 - 10/23 @ 2155 





WAVE PROPERTIES 


615 


10/23 (a). 2155 


a (1/m) 


.572645 


M. (m) 


- . 04948 


CT. (m) 


.275683 


CT (m) 


.267219 


d(m) 


1.74 


a/d 


. 153574 



Figure 4.1 EUstogram Data 
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Gage 635 - 10/23 @ 2155 



1.60 



1.40 



S 120 

e 

£ 1.00 
tfl 

1 0.80 
o 

= 0.60 
XI 

£ 

2 0.40 
0.20 



0.00 

o 
o» 

9 

Disptac«m«nt (m) 



j Histogram Data POF Gaussian 



Figure 4.2 Histogram Data 




Table 4.2 Gage 635 - 10/23 @ 2155 





WAVE PROPERTIES 


635 


10/23 ® 2155 


a (1/m) 


.481441 


H. (m) 


- .0547 


CT. (m) 


.318871 


a(m) 


.30989 


d(m) 


2.83 


a/d 


.109502 




Table 4.3 Gage 645 - 10/23 @ 2155 





WAVE PROPERTIES 


645 


10/23 ® 2155 


a (1/m) 


.200988 


P. (m) 


- . 03084 


<T. (m) 


.343446 


CT(m) 


.340708 


d(m) 


4 .31 


CT/d 


.079051 



Figure 4.3 Histogram Data 
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0«ge 666 - 10/23 e 2166 




Displacement (m) 



Hstog^amData PDF Gaussian 



Table 4.4 Gage 655 - 10/23 @ 2155 





WAVE PROPERTIES 


655 


10/23 (a). 2155 


a (1/m) 


.104437 


fu (m) 


- .01771 


CF. (m) 


.340056 


a(m) 


.339137 


d(m) 


6.0 


CT/d 


.056523 



Figure 4.4 Histogram Data 



Gage 666 - 10/23 @ 2166 




histogram Data PDF Gaussian 



Table 4.5 Gage 665 - 10/23 @ 2155 





WAVE PROPERTIES 


665 


10/23 (2) 2155 


a (1/m) 


.126266 


H. (m) 


- .01835 


CT. (m) 


.365155 


CT(m) 


.364236 


d(m) 


6.08 


CT/d 


.059907 



Figure 4.5 Histogram Data 
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12 



c 0.8 



Gage 676 - 10/23 @ 2166 



0.6 



5 0.4 
o 
£ 



0.2 




Qsplacement (m) 



Kstogram Data PDF Gaussian 



Figure 4.6 Histogram Data 



Table 4.6 Gage 675 - 10/23 @ 2155 





WAVE PROPERTIES 


675 


10/23 (2) 2155 


a (1/m) 


.094349 


(m) 


- .01014 


CT. (m) 


. 375221 


CT(m) 


.374948 


d(m) 


6.82 


a/d 


.054978 



Cage 626 • 10/23 @ 2166 




I E222ES3 Histogram Data PDF Gaussian 



Table 4.7 Gage 625 - 10/23 @ 2155 





WAVE PROPERTIES 


625 


10/23 ® 2155 


a (1/m) 


.069914 


M. (m) 


- .00819 


a. (m) 


.394571 


a(m) 


.394401 


d(m) 


9.66 


a/d 


. 040828 



Figure 4.7 Histogram Data 
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Table 4.8 Gage 615 - 10/25 @ 1315 





WAVE PROPERTIES 


615 


10/25 (2), 1315 


a (1/m) 


.372027 


(m) 


- .09863 


<T. (m) 


.498230 


a(m) 


.47981 


d(m) 


2.23 


a/d 


.215161 




Table 4.9 Gage 635 - 10/25 @ 1315 





WAVE PROPERTIES 


635 


10/25 (a). 1315 


a (1/m) 


.318011 


P. (m) 


‘ .10497 


CT, (m) 


.591069 


a (m) 


.573397 


d(m) 


3.32 


CT/d 


.17271 



Figure 4.9 Histogram Data 
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Figure 4.10 EUstogram Data 



Table 4.10 Gage 645-10/25 @ 1315 





WAVE PROPERTIES 


645 


10/25 (ffi 1315 


a (1/m) 


.256391 


ji. (m) 


- .10245 


CT. (m) 


.628243 


CT(m) 


.612186 


d(m) 


4.66 


CT/d 


.13137 




Table 4.11 Gage 655 -10/25 @ 1315 





WAVE PROPERTIES 


655 


10/25 (S>. 1315 


a (1/m) 


.269927 


P. (m) 


- .12356 


CT. (m) 


.658813 


a(m) 


.636821 


d (m) 


6.44 


a/d 


. 098885 



Figure 4.1 1 Histogram Data 
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Table 4.12 Gage 665 -10/25 @ 1315 





WAVE PROPERTIES 


665 


10/25 ® 1315 


a (1/m) 


.269743 


(m) 


- .12978 


or. (m) 


. 67434 


CT(in) 


,650698 


d(m) 


6.52 


a/d 


.09980 



Gage 676- 10/26 @ 1316 



0.7 




•2 - 1.2 - 0.4 0.4 1.2 2 2.8 3.6 



Displacement (m) 



Knissssa^ histogram Data PDF Gaussian 



Table 4.13 Gage 675 -10/25 @ 1315 





WAVE PROPERTIES 


675 


10/25 ® 1315 


a (l/m) 


.220143 


P. (m) 


- .12062 


CT. (m) 


.727669 


CT(m) 


.708477 


d(m) 


7.22 


cr/d 


.098127 



Figure 4.13 EUstogram Data 
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Table 4.14 Gage 625 -10/25 @ 1315 





WAVE PROPERTIES 


625 


10/25 (2) 1315 


a (1/m) 


.140009 


p. (m) 


- .10961 


CT. (m) 


.872244 


cT(in) 


.858792 


d(m) 


10.04 


ci/d 


.085537 



Based on a comparison with the water depth regions defined in Section 3.2. 1, it 
can be seen that the histogram data and the probability density fimction closely resemble 
the normal distribution even after the wave field enters intermediate water for milder sea 
conditions. This is shown in Figures 4.6 and 4.7. This impUes that the characteristics of 
the wave profile are dependent on the sea condition as well as the water depth. 

The histogram data for both the mild and severe sea conditions in Figures 4.1-4. 14 
show the distribution of the wave profiles progressively becomes more non-Gaussian as 
the waves propagate toward shallower depths. The non-Gaussian trend becomes more 
extreme as the sea condition becomes more severe. It is significant to note also that the 
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non-Gaussian characteristics are most pronounced at the shallower gage locations, Gages 
615 and 635. 



4.2 Gaiigsian/non-Gaussian Boundary Definition 

The histogram data shown in Figures 4. 1 - 4. 14 represent the same wave profiles 
depicted in Figures 3.2 - 3. 15. The transformation fi^om Gaussian to non-Gaussian 
profiles begin with subtle changes in the amplitude of the negative displacements and the 
width of the troughs. This initial transformation can be seen be comparing Figures 4.4 and 

t 

4.6. Figure 4.6 shows a nearly normal distribution. Figure 4.4 shows a slight elevation of 
the apex and a narrowing of the negative tail. It is also significant to note that the positive 
part of the distribution remains relatively unchanged. This indicates that the initial changes 
afiFect the trough region of the profiile more drastically than the positive region. Figure 4.4 
demonstrates that the trough region becomes shallower and wider. As the profile 
continues to shoal, changes begin to occur in the positive displacement region of the 
distribution. The peaks become narrower and higher. This can be seen graphically by 
comparing Figures 4.2 and 4.4. The positive region of the wave profile is not as full when 
compared with the normal distribution and the positive tail extends farther than the normal 
curve. 

As described above, significant transformations occur in the wave profile as the 
wave train propagates into shallower water. There is a point in this transformation when 
the wave profile deviates significantly enough fi'om the normal distribution that the waves 
can no longer be considered a Gaussian random process. Identification of this boundary is 
important to ensure assumptions made regarding the waves remain vaUd. Analysis of the 
histogram data shows that the non-Gaussian characteristics of coastal waves depends on 
water depth and sea severity. Based initially on analysis of the histogram data shown in 
Figures 4.5 and 4.6, it appears this transition boundary occurs when the non-dimensional. 
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wave standard deviation, a/d, is approximately 0.06. The non-dimensional parameter, 
a/d, is used because it considers both sea severity and water depth. Eight random sea 
conditions were selected having o/d values between 0.04 and 0. 10. These eight data sets 
are depicted in Figure 4. 15 in terms of ao as a function of a/d . The data points labeled A- 
H in Figure 4.23, each correspond to a distribution shown in Figures 4. 15-4.23. Data 
points A-G were intentionally grouped around the 10 meter water depth. This minimized 
any impact water depth and/or beach slope may have on the analysis results. The 
statistical data for these points are provided in Table 4. 15. 
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Figure 4.15 Sample Data Points to Analyze Gaussian/Non-Gaussian Boundary 
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Table 4.15 Statistical Properties of Data Points A-H of Figure 4.15 



POINT 


a/d 


H, 


Water Depth, d 


A 


0.1000 


3.08 m 


9.51 m 


B 


0.0800 


2.74 m 


9.67 m 


C 


0.0710 


2.47 m 


9.64 m 


D 


0.0616 


2.47 m 


10.03 m 


E 


0.0600 


2.43 m 


10.13 m 


F 


0.0520 


1.78 m 


8.90 m 


G 


0.0398 


1.57 m 


9.88 m 


H 


0.0590 


4.77 m 


20.20 m 



Comparison with Normal Distribution 
o/d = .100 




-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 

Displacement (m) 



Prob. Density Fen. Normal Distribution 



Figure 4.16 Point A of Figure 4.15 
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Comparison with Normal Distribution 
d/d = .080 




FVob. Density Fen. Normal Dstribution 



Figure 4.17 Point B of Figure 4.15 



Comparison with Normal Distribution 
aid = .071 




Rob. Density Fen. Nornrel Distributbn 



Figure 4.18 Point C of Figure 4.15 
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Comparison with Normal Distribution 

aid = .0616 




Rrob. Density Fen. Nornel Distribution 



Figure 4.19 Point D of Figure 4.15 




Figure 4.20 Point E of Figure 4.15 
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Comparison with Normal Distribution 
a/d = .052 




-2 - 1.5 -1 - 0.5 0 0.5 1 1.5 2 



Displacement (m) 



Rrob. Density Fen. Normal □stribution 



Figure 4.21 Point F of Figure 4.15 



Comparison with Normal Distribution 
a/d = .0398 




FVob. Density Fen. Normal Distributbn 



Figure 4.22 Point G of Figure 4.15 
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Comparison with Normal Distribution 
a/d = .059 




Displacement (m) 



Prob. Density Fen. Normal Dstribution 



Figure 4.23 Point H of Figure 4.15 



Figures 4. 16 - 4.23 confirm that the region between a/d = 0.06 and a/d = 0.08, 
shown in Figures 4.20 and 4. 17 respectively, appears to be when the wave profile 
transforms fi^om Gaussian to non-Gaussian. A closer inspection of this region is shown in 
Figures 4. 18 - 4.20. The distribution shows a considerable increase in deviation from the 
normal between a/d = 0.06 and a/d = 0.0616 as shown in Figures 4.20 and 4. 19 
respectively. 

Due to the subjective nature of this boundary, a general boundary can be defined 
based on probabihstic analysis. This boundary, defining the conditions under which the 
wave profile can no longer be considered Gaussian, occurs when a/d =0.06. 

A similar study was conducted by Ochi and Wang (1984); however, they defined 
the boundary condition between Gaussian and non-Gaussian wave profiles in terms of the 
wave record parameter, X^ represents the skewness of the wave profile which they 

determined to be the dominant parameter affecting the non-Gaussian characteristics of 
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coastal waves. Ochi and Wang concluded that coastal waves in seas for which X 3 is less 
than 0.2 can be considered as a Gaussian random process. For comparison purposes, this 
boundary condition is plotted along with the findings of this present study in Figure 4.24. 
In terms conq)atible with the approach of this present study, their boimdary condition can 
be defined as a /d = 0.0543. It is encouraging that there shows very close agreement 
between both findings. Additional analysis was contributed by Ochi and Ahn (1994b) 
whose data points are shown in Figure 4.24 as white circles. Ochi and Ahn (1994b) 
defined the transition criterion between Gaussian and non-Gaussian as occurring when a < 
0.01 where “a” is one of the parameters in the non-Gaussian probability density fimction. 
The white circle data points define the limiting cases based on his criterion. This third 
method for determining the boundary between Gaussian and non-Gaussian wave profiles 
agrees very well with the other two further verifying the vahdity of this finding. 
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Figure 4,24 Gaussian/Non-Gaussian Boundary Defined 



CHAPTERS 

TREND ANALYSIS OF PARAMETERS 
5.1 Relationship between o and Water Depth 

As was discussed earlier, the distribution of the wave profile is highly dependent 
on sea severity and water depth with sea severity being defined in terms of the wave 
profiles’ standard deviation, c. ’ Since these two terms will be used considerably in the 
following trend analyses, their relationship with each other will be examined. 

Figure 5. 1 shows examples of simultaneous measurements of o plotted as a 
fimction of distance offshore for two cases. It is evident fi^om Figure 5.1, that the 
standard deviation of the wave profile is influenced by water depth. As water depth 
decreases, the relative magnitude of the standard deviation decreases. It is also seen that 
for more severe sea conditions, represented by 10/25 @1315 data in Figure 5.1, the 
overall change in a as the wave field propagates into shallower waters is greater as 
compared with milder sea conditions (10/23 @ 2155). Based on theoretical 
considerations of a, it can be assumed that a approaches zero as distance oflFshore 
approaches zero. This trend is observed in Figure 5. 1 since both cases shown appear to 
converge to the origin. From the information on the surf zone locations for these two 
cases computed in Chapter 3, it appears that o begins to decrease rapidly towards zero 
once the break point is reached. 

The information shown in Figure 5. 1 based on two cases is shown in Figure 5.2 for 
all data sets analyzed in this study. Instead of a being plotted as a fimction of distance 
offshore. Figure 5.2 plots a as a function of water depth. This will take into consideration 
any local changes in mean sea level. 
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Plot of a VS Distance Offshore 




— 10/23 © 2155 10/25 © 1315 — Mean Water Depth Profile 



Figure 5.1 Spatial Comparison of a at Mild and Severe Sea Conditions 



The trends evident in Figure 5.2 reinforce the analysis done for the two cases in 
Figure 5. 1. The upper bounds of the discrete data represent the severe sea conditions 
while the lower bounds represents the mild sea conditions. 
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Figure 5.2 All a Data Plotted as a Function of Water Depth 
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It can also be inferred from Figure 5. 1 and Figure 5.2 that as the sea state becomes 
milder, the waves’ standard deviation approaches a constant value and hence the sea 
severity approaches a constant before decreasing towards the origin. 

In the absence of elaborate wave data, the value for a can be readily determined 

with a knowledge of the local and breaking wave heights and the relative water depth, 

defined by — ^ , where g = gravitational constant and T = wave period, using a method 

gT 

such as stream fimction theory (Dean 1974). To ensure the use of accurate wave height 
values in applying the stream fimction theory, nonlinear shoaling methods should be used 
in finite water depths to determine the breaking and local wave heights. 

5.2 Non-Gaussian PDF Parameter Analysis 

The non-Gaussian probability density fimction parameters, a, (i. and o,, will be 
analyzed in non-dimensional form spatially as a fimction of distance oflfehore, against each 
other in non-dimensional form, and with aa as a fimction of a/d. The first analysis of the 
non-Gaussian probabihty density fimction parameters will involve comparing the non- 
dimensional parameters against the distance the gage was located offshore. This will 
provide a spatial orientation of the parameters which will disregard any local fluctuations 
in water depth due to tidal cycle of storm surge. The results are shown in Figures 5. 3-5. 5. 

5.2.1 Spatial Analysis of Parameters 

AH three parameters shown in Figures 5. 3-5. 5 appear to gradually increase in 
absolute magnitude untU approximately 150 meters to 175 meters offshore. The mild 
condition represented by 10/23 @2155 consistently yields lower values when compared 
with the more severe sea condition on 10/25 @ 13 15. As the sea condition becomes more 
mild, the wave profile resembles the normal distribution. When the wave profile 
represents exactly the normal distribution, “a” equals zero and a» equals a. It can be seen 



Plot of Bottom Profile and "aa" Comparison 
Distance Offshore (m) 
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Figure 5.3 Spatial Comparison of aa at Mild and Severe Sea Conditions 
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Figure 5.4 Spatial Comparison of p./a at Mild and Severe Sea Conditions 
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Figure 5.5 Spatial Comparison of a./a at Mild and Severe Sea Conditions 
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from Figures 5.3 and 5.5 that as approaches zero under mild conditions and s./s 
approaches one. Figure 5.4 shows that m«/s approaches zero under these same conditions 
which is also consistent with the theoretical limit of the non-Gaussian probability density 
function. 

There is also a significant perturbation evident in Figures 5.3-5. 5 around 150 
meters ofishore. The mild sea condition shows this abrupt change occurring seaward 
when compared with its location imder severe sea conditions and the perturbation is more 
pronoimced in the severe sea condition case. In general, all trends analyzed for the non- 

t 

dimensional parameter, m»/s, ?re almost identical results to those found for as and s*/s but 
in the mirror image. This is due to the fact that the value of m. is negative. The negative 
sign is supported by the theory used to develop the probability density functions’ three 
governing equations. 

The next spatial analysis will consist of comparing the non-dimensional parameters 
with water depth. This comparison will be sensitive to any influence of local depth 
changes occurring from tidal cycles and/or storm surges. Figures 5.6 - 5.8 show the 
results of this analysis. 

The trends from this analysis are the same as those determined by Figures 5.3-5. 5. 
Due to the difference in scale of the abscissa the curves based on water depth appear to 
show much smoother transitions. Figures 5. 6-5. 8 include two additional cases 
representing intermediate sea conditions between those represented by 10/23 @2155 and 
10/25 @ 1315. These two additional cases approximate the trends displayed by the mild 
sea condition of 10/23 @2155. This is consistent with the data and time of these cases 
since they occurred prior to the severe storm so the sea condition was not yet impacted 
significantly. 

As was discussed earUer, in the deeper water conditions aU three parameters tend 
to show convergence towards their respective values for a Gaussian wave profile. In 
extremely shallow water conditions, the parameters also tend to converge towards a 
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Figure 5.6 ao as a Function of Water Depth 
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Figure 5.7 p./a as a Function of Water Depth 
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Figure 5.8 o./a as a Function of Water Depth 
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constant value \^dlich is different depending on the sea condition. For the milder sea 
conditions, the parameters tend to converge to the values shown in Table 5. 1. 



Table 5.1 General Shallow Water Limits of Parameters 



Non-Dimensional Parameter 


MUd - 10/23 (% 2155 


Severe - 10/25 (% 1315 


aa 


0.15 


0.18 


p*/a 


-0.18 


-0.22 to -0.25 


a,/a 


1.033 


1.04 to 1.045 



When the separate plots for each non-dimensional parameter are compared in both 
spatial cases, it is apparent that similar trend descriptions fit all three parameters. For this 
reason, a closer analysis will be made among the three variables to determine if this 
similarity can be quantified. The next trend analysis consists of comparing the three 
parameters with each other. 

5.2.2 Relationship Between the three non-dimensional parameters, aa. p,/a and a,/g 

As discussed earher, the three non-dimensional parameters appear to be highly 
correlated with each other. A closer look will be taken to see if specific relations can be 
developed connecting the three non-dimensional parameters together. Figure 5.9 and 
Figure 5. 10 show plots of the raw data representing p.,/a and a,/a compared with aa 
respectively at the extreme water depth locations, gage 615 and gage 625. 

All data sets from Gages 615 and 625 were plotted since they represent the 
shallowest and deepest gage locations. There are very definite relationships which hold 
constant between the extreme range of water depths and sea conditions, that link aa to 
p.»/a and a»/a. The frmctional relationship can be determined only after considering all 
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the data points from all water depth and sea severity conditions to ensure it remains 
consistently valid. A deterministic relationship between these non-dimensional parameters 
is a topic of active research. 



Trend of aa vs pJo 

0 0.05 0.1 0.15 0.2 

0 

-0.05 
- 0.1 

- 0.2 
-0.25 



# Gage 61 5 o Gage 625 

Figure 5.9 Plot of ixJc as a function of aa for Gages 615 and 625 
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Figure 5. 10 Plot of c./g as a function of ao for Gages 615 and 625 
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5.2.3 Relationship Between aa and a/d 

Based on the apparent correlation between the parameters of the non-Gaussian 
probability density function, aa will be used to represent the general trends of the other 
two parameters when plotted as a function of a/d. This relationship is provided in Figure 
5. 1 1 to demonstrate general trends between the wave aa, and the local environment 
represented in terms of sea severity/water depth in non-dimensional form. The lower 
values of a/d correspond to the deeper gage locations and the trend is for these values to 
approach the values defined by wave profiles having a normal distribution. This plot also 

t 

reinforces the analysis done in Chapter 4 which determined the boundary values of the 
Gaussian distribution. At a/d = 0.06, the aa values appear to begin converging towards 
the Gaussian limit located at the origin of Figure 5. 1 1. 

The higher values of a/d represent the shallower water depths. It appears that at 
these values aa tends to approach a finite asyn^totic limit. This is consistent with the 
findings fi'om the spatial analysis figures. In between these two limi ts of Figure 5.11 there 
is a depression in the severe sea condition case. This phenomenon must be analyzed more 
thoroughly to ascertain a reasonable explanation for its existence and will be the subject of 
a later study. 
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Figure 5.11 Plot of aa as a Function of o/d for Mild and Severe Sea Conditions 



CHAPTER 6 

SUMMARY AND CONCLUSIONS 



It is shown, based on the results of this study, that wind generated waves exhibit 
significant changes in their profile forms as they propagate into finite water depths. The 
actual nature of this transformation has been a subject that has not been covered in any 
significant detail. The unique wave data records collected during the ARSLOE Project of 
1980 enabled a detailed study on this topic to be produced. 

One objective of this study was to verify the appUcability of the non-Gaussian 
probability density fimction developed by Ochi and Ahn (1994b) defined as 

f(x) = -=L — exp — —(I - ya^. - exp(-rax)f - yax 

■yJln-cT, [ 2(yaa.) 

, \2^for.x>Q 
where, y = \ 

^ 3.00/orx<0 

Through statistical analysis of the ARSLOE data, it has been confirmed that this 
probabUity density function can be used with confidence to represent wave profiles in 
finite water depths. It correlates very well with histogram data over an extensive range of 
sea conditions. The distribution of the wave profile in finite water depths shows a 
dependence on water depth and sea severity. In particular, as water depth decreases, the 
distribution becomes increasingly non-Gaussian. Also, as sea severity increases the same 
trend occurs. 

With the appUcability of the probabiUty density function verified, an analysis of the 
boundary condition where the wave profile can no longer be considered Gaussian was 
conducted. A general relationship has been proven to exist between the largest significant 
wave height and the local water depth where the wave profile can still be considered 
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H O’ 

Gaussian. This relationship is defined as — — = 0.24 or — = 0.06. Both forms of this 



relationship are identical since the significant wave height, Hs is defined as four times the 
standard deviation, o, in deep water conditions. This relationship has been verified using 
three different approaches with each yielding similar results. In shallower water 
conditions, this boundary, defining the transition fi^om Gaussian to non-Gaussian wave 
profiles, appears to shows a dependence on beach slope. This is a qualitative observation 
fi-om the present study and additional research must be done to quantify this finding. This 
dependence on beach slope may be similar to the influence of beach slope on breaking 
waves. The linear result defined by this study is akin to the earliest breaking criterion 
based on Sohtary Wave Theory. 

At present, the only way this probabihty density fimction can be defined is through 
knowledge of statistical properties of the wave field of interest. In particular, the mean, 
variance and third moment of the wave profile must be known. This requires that wave 
profile data be obtained in order to evaluate the parameters of the probabUity density 
fimction. 

Another objective of this study was to determine relationships linking the 
parameters of the probabihty density fimction to local field conditions. If the wave profile 
statistics can be defined in terms of local environmental conditions such as sea state and 
water depth, the probabihty density fimction can be defined without wave data. This study 
has revealed quahtative trends between the parameters of the probabihty density fimction 
and local field conditions based on the extensive data comphed fi'om ARSLOE. From 
these quahtative trends, it is evident that the specific relationships between the parameters 
and environmental conditions are complex and may require consideration of additional 
parameters besides sea severity, a, and water depth. The pursuit of a fimctional 
relationship between the parameters and the local environment continues to be a topic of 
considerable interest and active research. 



APPENDIX 

COMPARISON OF ORIGINAL AND SMOOTHED “a” VALUES 



It can be seen that the nonlinearity parameter, “a” is characterized by a significant 
amount of local flucuations superimposed on the overall trend of the parmeter. Based on 
the definition of “a” as a nonlinearity parameter used to describe the probability density 

t 

fimction of a wave profile, these small flucuations can be smoothed without any noticable 
loss of accuracy. It is highly unlikely that the probability density fimction of a wave profile 
changes so sporadically over such short periods of time as indicated by the minor 
variations of “a”. There was not any one smoothing technique that was used in the 
smoothing of “a” at all gage positions. Instead, several smoothing techniques were 
evaluated and the technique which produced the lowest mean absolute deviation (MAD) 
was the technique used for that particular gage position. 
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Figure A. 2. Plol of Gage 635 Original and Sinootlied "a" Values 
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iMguie A. 4. Plot of Gage 655 Original and Sinootlied "a" Values 
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I'igme A. 5. Plot ol'Giige 665 Original and Smoothed "a" Values 
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I'igiire A. 6. Plot ofCjiigc 675 Original uiul SimxXlicd "a" Values 
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